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We investigate the creation of a relative phase between 
two Bose-Einstein condensates, initially in number states, by 
detection of atoms and show how the system approaches a co- 
herent state. Two very distinct time scales are found: one for 
the creation of the interference is of the order of the detection 
time for a few single atoms and another, for the preparation 
of coherent states, of the order of the detection time for a 
significant fraction of the total number of atoms. Approxi- 
mate analytic solutions are derived and compared with exact 
numerical results. 

PACS number(s): 03.75.Fi, 05.30. Jp, 42.50. Ar 



I. INTRODUCTION 

The first experimental realisations of Bose-Einstein 
condensation in dilute atomic gases have opened up 
a new field in atomic and quantum physics. Despite its 
apparent complexity, the condensate is well-described by 
a macroscopic wavefunction, or complex field, obeying a 
nonlinear wave equation (the Gross-Pitaevskii equation) . 
The nature of this complex field and in particular its ap- 
parent phase has been the subject of some debate. 

One specific topic which has been frequently investi- 
gated is the question of the coherence properties ||,|| of 
Bose-Einstein condensates as established in interference 
experiments |3]j7]]. The underlying question is whether 
the phase or, more precisely, the relative phase of two 
condensates is created by a spontaneously broken sym- 
metry H or by some other mechanisms It has been 
su ggest ed ]To| (and studied in detail by several authors 
|ll)-|l6||) that such a relative phase can be created by 
the detections of individual atoms in an interferometric 
setup, that is, where the origin of the detected atoms is 
intrinsically unknown. This process gives rise to a defi- 
nite (but unpredictible) relative phase of two condensates 
even if the initial states of the condensates are of unde- 
fined phase, as for initial number states fL0[[ll|jl3[ |, initial 
Poissonian states Jl^,0 , or initial thermal states H] , by 
entangling the states of the two condensates. 

In this work we will study in detail the creation of 
coherence between two condensates which initially have 
well defined occupation numbers, not only in the limit 
where the number of atoms detected is small compared 
to the total number of atoms but also in the long time 
limit. We concentrate thereby on two manifestations of 
coherence: the creation of interference fringes and the 
evolution of the compound two-condensate system to- 
wards a kind of coherent state. We start with an ideal- 



ized model (Sees. [II - IV) and later introduce more realis- 
tic features including atomic collisions (Sec. 0). We use 
exact numerical simulations combined with approximate 
analytical solutions to determine the time evolution of 
our system. We show that these two measures of coher- 
ence are created on very different time scales. The rela- 
tive phase associated with the appearance of interference 
fringes develops in a time of order l/(N r y), where N is 
the total initial number of atoms in the two condensates 
and 7 is the rate at which atoms leak out of the conden- 
sates and are detected. The approach towards coherent 
states, in contrast, takes a larger time of order I/7. These 
two timescales correspond to timescales recently identi- 
fied for the interaction of a Bose-Einstein condensate or 
other bosonic system with its environment [fl7|| . In this 
case, a well-specified atom number will change on the 
timescale l/(Nj) but any coherence present decays on a 
timescale I/7. 



II. MODEL 

Let us first introduce the model system which we use 
to investigate the coherence properties in Bose-Einstein 
condensation p~0} [l6[] . We consider two independent non- 
interacting single-mode Bose-Einstein condensates with 
creation (annihilation) operators (a) and (b), re- 
spectively. Atoms are leaking out of the condensates 
and are detected individually and spatially resolved. The 
same detectors simultaneously monitor decays from both 
condensates, that is, atoms coming from different con- 
densates are allowed to interfere, see Fig. 0. Thus, the 
detections are described by the annihilation operators 
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where <fi £ [—tt, tt] is related to the position of the detector 
x by <j) = px/fi where p is the momentum of the atoms 
leaking out of the condensates. 

Assume now that the system can be described at a cer- 
tain time t by a wave function \ip). Then the probability 
of detecting the next atom at position <j) is given by 



P{ct>)=N{TP\a{<t>ya{<f>m 



(2) 



where the normalisation constant J\f is chosen such that 
d<t>P(ct>) = l. (3) 



This probability function can be rewritten as 



1 



P(fa = — [l + c cos^-d)} 



where 
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(5) 



is the visibility of the interference fringes conditioned on 
the quantum state and 9 gives the most likely po- 
sition of detection of the next atom. It should be em- 
phasized that this visibility is not the one obtained by 
detecting a large number of atoms from the initial state 
IV'), but the one obtained by preparing this state very 
often and measuring a single atom in each run. This dif- 
ference is important since every detection changes the 
state of the system and thus changes the conditional vis- 
ibility P c . 
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FIG. 1. Schematic representation of the interfering 
Bose-Einstein condensates. 

This system has already been investigated by several 
authors JlOj— []~6|] in order to show that the detection of 
atoms breaks the underlying symmetry of the system and 
thus creates interference fringes. This is true even if the 
initial state of the system does not exhibit any prefered 
phase so that C = 0. It is the main purpose of our 
work to quantify this creation of coherence between two 
initially uncorrelated Bose-Einstein condensates. 



III. CREATING INTERFERENCE FROM 
INITIAL NUMBER STATES 

In this section we will discuss the creation of interfer- 
ence fringes as a consequence of consecutive detections of 
atoms when the two condensates are initially in number 
states. The full system is given initially by the quantum 
state 



ni,n 2 ) 



(6) 



After the detection of k atoms at positions fa, fa, fa 
the (unnormalizcd) state of the system is then 



|^ fc ) = (a + 6e^)...(a + 6e^ 1 )|Vo). 



(7) 



The conditional probability of detecting the (k + l)th 
atom at fa-+i then reads 

P(fa+i\fa, -,fa) 



2tt (ip k \tfa + tfb\ip k ) 

1 (fa +1 \ip k+ i) 
2ir (N - k){ip k \ip k ) 



(8) 



where N = n\ + n 2 - Thus, the probability for the se- 
quence of detections fa, fa, fa is 



P(fa,...,fa) = P(fa)P(fa\fa). 

_ i (fa\i>k) 



(2Tr) k N(N — 1)...(N -k + 1)' 



(9) 



The conditional visibility (3 (we will write instead of (3 C 
in this section in order to simplify the notation) for the 
state after k detections is 



= 



2\(i, k \a)b\^ k )\ _ 2\(fa\a?b\4> k )\ 



(fa\rfa + tfb\fa) (N-k)(ip k \i/j k )' 
Thus the average visibility after k detections is 



(10) 
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Note, however, that this mean conditional visibility is 
rather difficult to access experimentally since it involves 
the averaging over ensembles of experiments where each 
ensemble consists of repeatedly preparing the quantum 
state of the system by the same sequence of detections 
from the same initial number state and measuring the 
position of the next detection. Nevertheless, this proves 
to be a useful measure theoretically especially to describe 
the time scale on which interference is created as we will 
show later in this section. 

For the first two detections, that is, k = 1, 2, the 
integral in Eq. ( |TT| ) can be evaluated analytically with 
the results 



(0) 



2n\ni 



{ni + n 2 ) 2 - (ni + n 2 ) ' 

(0h = -{0)i. 
ir 

For m =n 2 = N/2 we thus obtain 

OTi = 1 1 



2 1 - 1/N 



(12) 
(13) 

(14) 



For any initial number of atoms this is larger than 1/2, 
meaning that the first detection already increases the av- 
erage visibility from zero to more than half its maximum 
value of one. 

For k > 2 Eq. ( |TT| ) can no longer be evaluated analyt- 
ically but we will derive an approximate solution in the 
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following. Let us assume that the system starts in the 
quantum state \4>o) = \n,n) (same number of atoms in 
both condensates) and all detections occur at the same 
position <p. Of course, this is a highly unlikely detection 
sequence, but this assumption gives surprisingly good ap- 
proximate results as we will see. Without loss of gener- 
ality we may assume <fr = 0. Thus 

|V^)-(a + 6) fe |Vo). (15) 
From this we obtain the norm of the state as 
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We now approximate the binomial coefficients by Gaus- 
sians using 



-f(m-fe/2) 2 



(17) 



and replace the sum over m by an integral from — oo to 
oo which yields 
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(2n - fc)! 7T ^k{An - fc) 
Analogously we obtain 



(18) 
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(2n-fc-l)! 71" v /fc(4n - k - 2) 

(19) 

and thus from Eq. (|J) the final result 

(3 k « e- 1 / fc (20) 

independent of n (to order 1/n). Hence, independent 
of the initial number of atoms in the condensates, it al- 
ways needs the same (and very small) number of detected 
atoms to create the interference. If the initial total num- 
ber of atoms is N and each atom decays with a rate 7 
out of the condensate, then the time to create the inter- 
ference pattern will thus be of the order of l/(N r y). 

In Fig. ||we compare the approximate result of Eq. (|2(i| ) 
with the exact numerical solution of Eq. ([To]) . The com- 
parison shows that the approximation is excellent for k 
as low as 2 (10% difference), 3 (5%) and 4 (3%) for 
n = 100, even if the approximation made in Eq. Jl7j ) 
only holds for k 3> 1 . We also plot the numerical solution 
of Eq. (|l|), that is, the visibility averaged over all pos- 
sible outcomes with the appropriate probabilities which 
we obtained by Monte-Carlo simulations of the detection 
process. We see that the difference in the creation of in- 
terference fringes between the case where all atoms are 



detected at the same place and the case where all pos- 
sible detection positions are allowed is relatively small. 
This is somewhat surprising given the great difference 
between the extremely peaked distribution of detection 
positions of our approximation compared with the ex- 
pected sinusoidal behavior [see Eq. (|)]. It is, however, 
a consequence of the fact that only a small number of 
detections determine the interference pattern for all sub- 
sequent detections. 
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FIG. 2. Average conditional visibility (/3) vs number of de- 
tected atoms fc: exact numerical solution (solid curve), nu- 
merical solution if all atoms are detected at the same posi- 
tion (dashed curve), approximate analytical solution (dotted 
curve). The initial state is m = na = 100. 

So far we have restricted ourselves to the case of equal 
initial occupation numbers of the two condensates. In 
the case of unequal initial numbers the approximation 
assuming that all the detections occur at the same po- 
sition fails as it turns out that this changes the relative 
occupation of the two condensates (which is constant if 
the condensate decay rates are equal). An analytic ap- 
proximation can be found, however, assuming that the 
number of detected atoms k is small compared to the to- 
tal number of atoms N. In this case the normalisation of 
the state after k equal detections is 



(tpk\ipk) = 



k 2 

E(m) (ni,n 2 |(a t ) ro a m (6 t ) fc - m 6 fc - m |n 1 ,n a ) 



m=0 
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(21) 



The latter expression can be approximated by replacing 
the binomials with Gaussians and the sum by an integral 
as before. Applying the same procedure to (iph\a'b\ipk) 
one finally obtains the conditional visibility 



n\ + ri2 



(22) 



the long-time limit of which (fc ^ 1) has also been found 
in Ref. @. 
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Three features are worth mentioning here. First, the 
evolution of the visibility as a function of the number of 
detected atoms is the same as in the case of equal initial 
atom numbers. Hence, also in this case only a few detec- 
tions are required to establish the interference pattern 
and hence the corresponding time scale is again given 
by l/(Nj). Second, the maximum possible visibility is 
reduced to a value which depends on the ratio of the ini- 
tial occupation numbers of the two condensates. Third, 
this maximum visibility is exactly the same as for initial 
coherent states of the condensates with the same mean 
numbers of atoms. 
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FIG. 3. Average conditional visibility (/3) vs number ri2 of 
atoms initially in one of the condensates, the total number 
is ni + 7i2 = 100. Solid curve: exact numerical result after 
99 detections, dashed: exact result after 5 detections, dotted: 
approximate analytic result after 99 detections, dash-dotted: 
approximate result after 5 detections. 

In Fig. H we compare these results with the exact nu- 
merical solutions for arbitrary detection positions. We 
see that in the long-time limit, with nearly all of the 
atoms detected, the agreement is exact, thus the visibil- 
ity approaches the value found for coherent states of the 
same mean atom number. After a small number k of 
detections a more significant deviation from our approx- 
imate result is found, especially for significantly differing 
initial occupation numbers ri\ and ri2- However, we still 
find that Eq. ( p2] ) predicts the correct order of magnitude 
for the time required to establish the interference. 



IV. PREPARATION OF COHERENT STATES 
FROM INITIAL NUMBER STATES 

In this section we will discuss how the system state 
approaches a coherent state in course of the sequence of 
detections. However, in order to simplify the discussion 
and the numerical simulations we will assume that the 
total number of atoms in the two condensates is exactly 
known at any time; the system is initially in a number 
state |m, ri2) with ri\ atoms in condensate a and ri2 atoms 
in condensate b and the number k of detected atoms is 
known at any time. Thus the system will never approach 



a harmonic oscillator coherent state, that is, an eigen- 
state to the annihilation operators a and b, since these 
are superposition states of different total atom numbers. 
Instead it will be shown that the system approaches a 
state which can be described as the restriction of a co- 
herent state to the subset of states with a fixed total 
number of atoms. Hereafter we will refer to these states 
as "atomic coherent states" as they are a representation 
of the atomic coherent states . We can define these 
states to be 
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(a*n + tfv) N \0,0) 



(23) 



where /i and v obey the relation |/i| 2 + |i/| 2 = 1. Some 
properties of these atomic coherent states and their rela- 
tion to the coherent states are discussed in Appendix [a|. 

The specific case of these states with /j = e %< ^ j\p2 and 
v = e~ 1 ^ /V2 has also been refered to as "phase states" 
PPL 
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( a t e # + b^e-^) N \0,0) (24) 
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since their conditional visibility is (3 C = 1. 

The problem which we will consider in this section is 
the following. Assume the system is initially in the num- 
ber state |^o) = 1^1)^2)- Then k atoms are detected at 
positions 0i, 02, <Pk and the resulting state \ipk) is 
analyzed. What is the probability of finding an atomic 
coherent state v)N-k, where N = ni+n2? To answer 
this we have to evaluate the probability function 



or equivalcntly 



K^feJ/X, v) N . 
(ipk\ipk) 



\(lp k \<P)N-k\ 
(ipk\ipk) 



(25) 



(26) 



in the case of rii = — N/2. 

Let us consider first the case of equal initial atom num- 
bers ni = ri2 — n and assume without loss of generality 
that the first atom is detected at position <pi = so that 
\tpx) — ( a + b)\n, n). We then find 



2n 
n 



(1 + cos 20) 



1 



(1 + cos 20) 



(27) 



where for the last approximation we have again used 
Eq. (0). One can easily check by numerical simulation 
that this approximation is highly acurate even for just 
a few atoms in the condensates. Hence the overlap of 
the state after one detection with any phase state is very 
small of the order of 1/y/n. 

For the analysis of the state after k detections we will 
assume, as in the previous section, that all atoms are de- 
tected in the same position, so that \tpk) = (a + b) k \ipo), 
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and will compare the results for P(4>) with numerical sim- 
ulations at the end. The state overlap after k detections 
is then 



(ipk\4>)2n-k 
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v /2 2 "- fe (2n- jfc)! 
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(28) 



which, after applying the approximation of Eq. (|l7|), be- 
comes 



{^k\4>) 2n- 
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n\2 2n ^ H2 n -k) 

e~ — — — . (29) 



This, together with the normalization factor, Eq. (18) 
yields 



„, ,, 1 k(4n—k) 4> 2 fctgn-fc) 

P(4>)~-\1— — n — -e 2 



(30) 



In the limit of k -C n this result is consistent with the 
results presented in Refs. [Ojiajla]. 
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FIG. 4. Probability P(0) of finding the phase state \ tj>) after 
fc atom detections from initial state |m = 100, n<2 = 100). The 
solid curves correspond to exact numerical solutions (averaged 
over 1000 Monte-Carlo simulations) and the dashed curves to 
the analytic approximation given in Eq. ({jo|). 

In Fig. |] we compare this analytic approximation for 
the probability function P(<f>) with the exact numerical 
solution obtained by Monte-Carlo simulations where all 
possible detection positions are taken into account. (On 
the scale of Fig. || the curves for the approximate solu- 
tion and the exact solution if all atoms are detected at the 
same position coincide almost exactly.) After only a few 
detections (fc = 10 in the figure) the probability function 
is already well approximated by a broad Gaussian with a 
relatively small maximum, so the overlap with any phase 
state is still small at this time. However, we note that 



the maximum overlap is larger than predicted by the ap- 
proximate analytic solution, which was derived under the 
assumption that all atoms are detected at the same posi- 
tion. This may seem surprising but can be explained by 
the fact that such a highly peaked position distribution 
of the detected atoms is far from the one expected for a 
coherent state. After the detection of half of the atoms 
(k = 100) the maximum overlap with a phase state is 
already close to one and the width of the Gaussian has 
decreased significantly. For a larger number of detected 
atoms (fc = 190) the maximum overlap still increases 
but the width of the probability function increases again 
which is due to the changing non-orthogonality o f th e 
phase states for changing number of atoms, see Eq. 
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FIG. 5. Maximum of the probability function P(/i, v) 
after k atom detections out of the N = rij + ni ini- 
tial atoms. The solid curves correspond to the initial 
state |-i/>o) = \n\ = 100, 712 = 100), the dashed curve to 
l^o ) — l n i = 200, ri2 = 50), and the dotted curve is the maxi- 
mum of the approximate analytic solution, Eq. (|o|). The first 
two cases are the results obtained from averaging over 2000 
Monte-Carlo simulations. 

In Fig. ^| we plot the maximum of the probability func- 
tion P(/i, v) as a function of the fraction of the detected 
atoms. The approximate solution given by Eq. (|3^) is in- 
dependent of the total number N = n% +ri2 of atoms ini- 
tially in the two condensates and is a quarter of a circle as 
a function of k/N. As already seen above, the approach 
to an atomic coherent state starting from a pure number 
state is in fact faster than given by the approximation. 
Here we also plotted the numerical result for an initial 
state with unequal atom numbers in the two condensates 
and we note that also in this case the approximation by 
Eq. ( pfj| ) is a relatively good one. 

Hence, we find that a significant fraction of the total 
number of atoms must be detected in order that the sys- 
tem approaches an atomic coherent state. Thus the time- 
scale for the preparation of such a state by detections is 
given by I/7 which widely differs from the time scale of 
l/(N"f) found in the previous section to be relevant for 
the creation of interference fringes. 
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V. GENERALISATIONS OF THE MODEL 



B. Proper time evolution 



A. Imperfect detection 

We will now generalize our results from the previous 
sections to the case of imperfect atom detection. Assum- 
ing that the detector efficiency is rj < 1 we expect that 
after k atoms have been lost from the condensates only rjk 
out of these are detected in the interferometric setup and 
thus contribute to the build-up of coherence between the 
two condensates, whereas the remaining (1 — rf)k atoms 
are simple losses from either of the two condensates. 

Hence, under the assumption that all detected atoms 
are found at the same position (as in the previous sec- 
tions) and the undetected atoms are coming with the 
same probability from either of the two condensates, the 
state after k atoms have been lost from the condensates 
is approximately 

m = (a + by k a^ k \n,n) 



(n-£k)\ 



(a + by"'\n-£k,n-£k) (31) 



[where £ = (1 — rj)/2] instead of the one given in Eq. (15). 
Thus, in our earlier results, Eqs. ( |20| ) and (|3C|), we only 
have to substitute n — > n — £fc and k — > rjk to obtain the 
approximate results for imperfect detector efficiency 



0k 



P(<j>) 



-i/( v k) 



(32) 
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rjk 



In — k + rjk 



t> fe=ERt . (33) 



As in the previous sections, comparison of these analytic 
approximations with exact numerical solutions shows ex- 
cellent agreement for the visibility and an actually 
faster approach towards coherent states as predicted by 
this approximation for P((j>). 

Eq. (^) shows that the only effect of imperfect de- 
tection on the creation of interference fringes is that the 
number of atoms detected per unit time is decreased and 
and thus it takes more time to detect the same number 
of atoms as for perfect detectors. The effect of losses of 
atoms from individual condensates does not seem to have 
any influence even if this changes the relative atom num- 
ber in the two condensates. However, since only a few 
atoms need to be detected in order to build up the in- 
terference, the fluctuations of the condensate occupation 
numbers remain very small compared to the total number 
of atoms and thus the maximum visibility is very close 
to the one for the unperturbed initial state. 

The situation is more complicated for the approach of 
the system state towards an atomic coherent state (or a 
phase state) because, not only the number of detected 
atoms, but also the total number of atoms left in the 
system plays an important role, so that P((j>) depends on 
k and n. The preparation of an atomic coherent state 
also occurs on a larger time scale than for r\ = 1 but the 
dependence on r\ is not simple. 



Another possible generalisation of our model is to in- 
vestigate the system evolution as a function of time in- 
stead of the number of detected atoms. Given a constant 
decay rate of 7 for individual atoms from the conden- 
sates the actual decay processes still occur in a proba- 
bilistic manner. Thus after a certain amount of time the 
number of atoms remaining in the system is uncertain. 
However, since the condensate decay follows an exponen- 
tial law the mean number of remaining atoms is 

(n)(t) - Ae" 27 * (34) 

and so the mean number of detected atoms is 



(k)(t) = N(l-e- 2 ^) 



(35) 



We can generalise the results of the previous sections to 
incorporate this time evolution by simply replacing the 
number k of detected atoms in Eqs. ( p0| ) and ([30|) by its 
mean value according to Eq. (|35|). 

Using this assumption we compared the exact numer- 
ical results for the cases of the system evolution versus 
number of detections and versus time. We found that the 
visibility is established slightly more slowly in the latter 
case. Let us consider, for example, the state of the sys- 
tem after such a time t that (k)(t) = 1. At this time 
there is still a significant probability that no atom was 
detected which greatly decreases the average visibility. 
On the other hand there is a certain probability that two 
or three atoms were detected which increases the average 
visibility, but since the difference between 0q and 0\ is 
much larger than the difference between fi\ and 02 [see 
Eqs. ([l2|), ([j~3|)1 the former term dominates and the aver- 
age visibility is smaller than /3i. However, this difference 
between the results for the two models decreases when 
more atoms are detected and, starting from 100 atoms 
in each of the condensates, after k — 10 detections the 
numerically found difference is already down to 1.7%. A 
similar agreement is also found for the maximum overlap 
of the system state with atomic coherent states in the 
two models. We may thus conclude that our approxi- 
mate analytic solutions also describe the time evolution 
of the system if one substitutes k — > (k)(t). 



C. Effect of collisions 

So far we have assumed the idealized case of noninter- 
acting particles in the condensates, so that our system 
was completely analogous to a system of photons in two 
high-quality cavities from which the photons decay and 
interfere. All of our previous results apply to this case as 
well || . 

However, it is well known that atomic collisions play 
an important role in the context of Bose-Einstein con- 
densation |I 20 2l[] . We will discuss the modifications to 
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our results in this case in the following. To this end, the 
free time evolution of the system between two quantum 
jumps now has to be replaced by the time evolution due 
to the collisional Hamiltonian which in our simple model 
of two single-mode condensates is njl9,p2] 



H = k [(a^a) 2 + (tfbf 



(36) 



The action of this Hamiltonian is to give different time 
dependent phases to the various number states of the 
quantum state of the system of, for example, an atomic 
coherent state. This dephasing gives rise to a time depen- 
dent "decay" of the coherence and therefore of the condi- 
tional visibility (3. This decay of the coherence counter- 
acts the creation of coherence due to atom detections and 
thus prevents the system of reaching a state of maximum 
visibility. 




FIG. 6. Time evolution of the conditional visibility /3 for 
different quantum states including collisions of atoms in the 
condensates. The initial quantum state is a state after 50 
detections from the number state [100, 100) (solid line), the 
state after 50 equal detections (a + £>) 50 1 100, 100) (dashed), 
and a phase state \4>)n=i50 (dotted), respectively. 

First we study the evolution of various quantum states 
under the action of the Hamiltonian ( |3"6| ) without atomic 
decays. Let the system initially be in an atomic coherent 
state with equal mean atom number N in the two conden- 
sates, IV') = \fJ-)fJ>)_N- The time evolution of the visibility 
(3, given by Eq. (H) can then be evaluated analytically as 



P(t) = [cos(2/t*)] 



N-l 



3 -2K 2 t 2 (Af-l) 



(37) 



where the latter approximation holds for small times t -C 
1/k. Note that the exact result of Eq. ( |37|) predicts the 
well known [|l3 16 2^,|23| revivals of the visibility after 
times which are multiples of 7t/(2k). 

As noted, however, in the preceding sections, the 
atomic coherent states are in general not a good approxi- 
mation to the state of the system until a significant num- 
ber of atoms have been detected. A better approxima- 
tion is provided by considering the evolution of the initial 



state \4>k) = (a 
tion of Eq. (FT' 



Using again the approxima- 



we obtain for this case 



e 1/fe exp 



-2nH 



2.2 k(2n 



k-1) 



An-k-2 



(38) 



In the limit of k <C n, we thus find that (3 decays as 
cxp(— K 2 t 2 k) which is much slower than the decay for an 
atomic coherent state (^). For k — > 2n the two expres- 
sions converge since the state \ipk) approaches an atomic 
coherent state in this case. 

In Fig. H we compare (3{t) for an atomic coherent state, 
the state \ipk), and the numerical result for a state after 
50 detections from an initial number state 1 100, 100) (all 
of these contain a total number of 150 atoms). We see 
that in any case the collision-induced decay can be well 
described by a Gaussian. The decay obtained for the 
state with simulated detections is faster than the one of 
( |38| ) which agrees with our finding of Sec. HI that 



Eq. 



an atomic coherent state is approached faster with arbi- 
trary detections than if all detections occur at the same 
position. 

We will now use these results to derive an approximate 
analytic expression for the visibility (3 after k detections 
including the effects of atomic collisions. Let us assume 
that at a given time t exactly k atoms have been de- 
tected from an initial number state \n,n} and that the 
visibility is /3 = {3q. Then the probability of detecting 
the next atom at time to + t is given by 



P(t) = 2n o7 e- 2no7t 



(39) 



where uq = 2n — k is the number of atoms left in the 
system at time to. Thus, if we write the decay of the 
visibility as 



0{to+t) = 0(to)e 



-t 2 /r 2 



(40) 



where r is given in Eq. (|3§|), then the visibility at the 
time immediately before the next atom detection is on 
average given by 

poo 

(*+()) = / dt (3(t + t)P{t) 
Jo 

= /3(i )2n o7 T^e"^ 2T2 [1 - $(n 0l r)} , (41) 

where $ denotes the error function. Let us now assume 
that the system is in steady state between the creation 
and the decay of (3, that is, the following detection in- 
creases the visibility again to its value (3(to) at time t . 
Hence, writing 



! 1 

fc 



(42) 



and using our previous result for the increase of (3 with 
the number of detections, Eq. (^p|), we obtain the follow- 
ing condition 



l-fWto+*)> = 1-TT^T (43) 
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with the solution 

k = 1 + 



1 



y/l-(P(t +t)) 

and thus the steady state visibility is 

1 



A- 



l + y/l-(P(to+t))' 



(44) 



(45) 




40 80 120 160 200 
k 

FIG. 7. Visibility f3 vs number of detected atoms k from 
an initial number state 1 100, 100) including atom collisions. 
Solid lines are exact results (averages of 1000 Monte-Carlo 
simulations), dashed lines are the corresponding analytic ap- 
proximations given by Eq. (fl5|). The parameters are (from 
top to bottom) k = 0.5 7, ft = 27, and k = 57. 

We compare this result with the results of Monte-Carlo 
simulations in Fig. ^ for different values of n. From this 
we note that the approximation yields values of that 
are too large. This is especially true for small numbers 
k of detected atoms (k < n). This was to be expected 
since (i) the approximation was based on the assumption 
of a steady state whereas it takes a certain number of de- 
tections to reach this state in the simulations and (ii) in 
the discussion of Fig. ^ we have already seen that the ac- 
tual decay of the visibility occurs faster than predicted by 
Eq. (H|) which in turn decreases the steady state value. 
On the other hand, the approximation yields values of (3 
below the numerical results for values of j3 < 0.8 where 
we know that the increase of the visibility due to detec- 
tions is faster than given by Eq. (p0|). 

Finally we will investigate the effects of atomic colli- 
sions in the condensates on the creation of atomic co- 
herent states. To this end we calculate the maximum 
overlap of the time evolved state 



\Mt)) 



-iHt 



\rl>k) 



(46) 



where \ipk) is the state after k equal detections from the 
number state |n, n) as used before, with the atomic co- 
herent states. Following the same steps as in Sec. IV one 
obtains 



max{P(0,f)} 



max{P((M = 0)} 
1 + [tnk(2n - k)/{2n)f 



(47) 



where max{P(</>, t — 0)} is the maximum overlap at time 
t = given by Eq. (|3p|). Hence, for k n this maximum 
overlap decays on a time scale of t ~ l/(nk) which is 
much faster than the time scale of the decay of the visi- 
bility where t ~ l/(nVk), Eq. ([||. This, together with 
the much longer time scale necessary to create a signifi- 
cant overlap with coherent states, shows that atomic col- 
lisions have a much larger effect on P{4>, t) than on the 
visibility. The time scales for the decay and the build-up 
of this state overlap also prevent the system from reach- 
ing a steady state and thus no analytic approximation 
analogous to Eq. (^5|) can be found. Hence we have to 
rely on numerical simulations in this case, see Fig. B. 
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FIG. 8. Maximum overlap max{P(/i, v)} vs number of de- 
tected atoms k from an initial number state |m = ni = 100) 
including atom collisions with collision rate k = 0.1 7 (solid 
line), ft = O.57 (dashed), and k = 57 (dotted). Each curve is 
obtained from averaging over 1000 Monte-Carlo simulations. 

We note that for small numbers of detected atoms the 
maximum overlap increases close to the case of n — 0. For 
larger k, corresponding to a broader distribution of the 
relative atom number between the two condensates, the 
dephasing due to the atom collisions starts to dominate 
and significantly reduces max{P(/i, v)} even for small 
values of k. However, when most of the atoms have been 
detected (k — > 2n) the function approaches unity as all 
one-particle states and the final vacuum state are exact 
atomic coherent states. 



VI. CONCLUSIONS 

In this work we have studied in detail the creation 
of coherence between two initially uncorrelated Bose- 
Einstein condensates according to the detections of in- 
dividual atoms in an interferometric setup. Our main 
finding is that first order coherence, as observed in the 
interference fringes of the condensates, is established on 
a time scale corresponding to the detection of only a few 
atoms, that is, within times t ~ l/(Nj), where TV is the 
total number of atoms in the condensates and 7 is the 
single-atom decay rate. On the other hand, high-order 
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coherence, which in this article we describe by the over- 
lap of the quantum state of the system with coherent 
states, is created by detecting a certain (large) fraction 
of the condensed atoms, so the time scale for this de- 
velop is given by t ~ I/7. These results were obtained 
from approximate analytic solutions as well as from exact 
numerical simulations. 

We have also investigated the effect of atomic collisions 
on these results. In this case the dephasing of the quan- 
tum state of the system due to the collisions counteracts 
the entangling effect of the atom detections. However, ac- 
cording to the widely different time scales a good visibil- 
ity of the interference fringes can be maintained whereas 
the atomic collisions effectively prevent the preparation 
of a coherent state even for relatively small collision rates. 
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APPENDIX A: ATOMIC COHERENT STATES 

The atomic coherent states Jl8| were defined in 
Eq. <M) to be 



f 



(Al) 



v'|2 



(A5) 
(A6) 



Clearly, the restriction of a two-mode coherent state to 
the N particle subset of states is an atomic coherent 
state. 

The atomic coherent states satisfy the following useful 
identities 

(A7) 
(A8) 

From these it immediately follows that the expectation 
values of a and b in an atomic coherent state are zero so 
that neither condensate exhibits a prefered phase. It is 
also clear that the mean number of atoms in the conden- 
sates a and b are N\/i\ 2 and iV|i/| 2 , respectively, and that 
the expectation value of a^b is N It follows that the 
conditional visibility for an atomic coherent state is 



a\fi,u) N = vJV/i|/z,^)jv-i, 
b\n,v) N = \/Nv\ii,v)n-\. 



P c = 2\ti\W\ 



(A9) 



which assumes its maximum value of unity if and only 
if the two condensates have the same mean occupation 
number. Restricting considerations to equal mean occu- 
pation numbers leads to the phase states |l^,[li| defined 
in Eq. (24). Further properties of the atomic coherent 
states can be found in the literature (l8|. 



where |/i| 2 + |^| 2 = 1. These are states with precisely N 
atoms shared between the two condensates and can be 
expressed as an entangled superposition of all product 
number states for which the total number is N: 

\^) N = ^f2( I l) 1/ \^y\n,N-n). (A2) 

The number of atoms in one of the condensates, there- 
fore, has a binomial distribution. Different atomic coher- 
ent states to the same number of atoms N are in general 
not orthogonal to each other, but have a nonvanishing 
scalar product 



N {[j,', v'lfj,, v) N = (fifj,'* + vv'*) N . 



(A3) 



The atomic coherent states are related to the familiar 
two-mode coherent states \a, a'} by 



\a,a') 



-(|a| 2 + |a'| 2 )/2 



E 

JV=0 



(14 



a- 



t\2\ N 



N\ 



(A4) 
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